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System and Method for 
Solving a Large System of Dense Linear Equations 

BACKGROUND OF THE INVENTION EL977166451US 

1. Technical Field 

5 The present invention relates in general to a system 

and method for solving a large system of dense linear 
equations. In particular, the present invention relates to 
a system and a method for solving a large system of linear 
equations using a system having multiple processing units 
10 with a common memory for sharing data. 

2. Description of the Related Art 

Many problems from many disciplines can be modeled 
using a system of linear equations. Linear equations, for 
example, may be used to obtain solutions for problems in 
15 physics, chemistry, engineering, computer science, etc. 

Small systems of linear equations (systems with only a 
few variables) can be solved by eliminating all but one 
variable, for which a solution is obtained. Variable 
elimination involves multiplying each equation by a 
20 constant or adding or subtracting one equation from 

another. Once one variable has been determined, the other 
unknown variables may be computed by back substitution. 

A system of linear equations has a corresponding 
matrix equation. The solution to this equivalent matrix 
2 5 equation is the same as the solution to the system of 

linear equations, and therefore, a solution to a system of 
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linear equations may be obtained by obtaining a solution to 
the equivalent matrix equation. 

A large system of linear equations is typically solved 
by obtaining a solution to the equivalent matrix equation. 
5 Numerically, it is more efficient to obtain a solution to 

the matrix equation. Efficient methods for solving the 
matrix equation include, for example, LU decomposition and 
Gauss elimination. The processes involve the exchanging of 
large amount of data between the processors since the 
10 results are interdependent. 

Thus, there is a need for a fast and accurate system 
and method for solving large systems of dense linear 
equations using multiple processors. The system and method 
should provide for fast and easy sharing of data results 
15 between the different processors to avoid the inefficient 

transmissions of large amount of data. 



Docket No. AUS920030703US1 



Atty. Ref. No. IBM-3064 



SUMMARY 

It has been discovered that the aforementioned 
challenges can be addressed by a method and a system for 
solving a system of linear equations using a system having 
a processing unit and one or more secondary processing 
units with a common memory for sharing data. The secondary 
processors share the processing involved 

A set of coefficients corresponding to a system of 
linear equations is received. The solution to the system 
of linear equations is the same as the solution to a matrix 
equation where the nxn matrix is formed by the coefficients 
of the unknown variables of the linear equations. To solve 
the matrix equation, an (n + l)xn augmented matrix may be 
formed by combining the nxn matrix with a vector of 
constant terms in the system of linear equations. 

The matrix is then divided into a plurality of 
coefficient blocks (sub-matrices), which are loaded into a 
common memory shared by a plurality of processors. Each of 
the processors may be programmed to apply matrix 
transformations to the coefficient blocks and return the 
sub-results to the common memory. 

A table containing a list of the matrix operations may 
be created in the common memory in order to monitor which 
operations have been performed and which are still pending. 
The table may be used to assign blocks requiring operations 
applied to available processors. The processors may 
indicate, using the table, the beginning and ending of a 
set of matrix transformations to a block of coefficients. 
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The foregoing is a suiranary and thus contains, by 
necessity, simplifications, generalizations, and omissions 
of detail; consequently, those skilled in the art will 
appreciate that the summary is illustrative only and is not 
intended to be in any way limiting. Other aspects, 
inventive features, and advantages of the present 
invention, as defined solely by the claims, will become 
apparent in the non-limiting detailed description set forth 
below. 



10 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention may be better understood, and 
its numerous objects, features, and advantages made 
apparent to those skilled in the art by referencing the 
accompanying drawings. The use of the same reference 
symbols in different drawings indicates similar or 
identical items. 

Figure 1 illustrates -the overall architecture of a 
computer network in accordance with the present invention; 

Figure 2 is a diagram illustrating the structure of a 
processing unit (PU) in accordance with the present 
invention; 

Figure 3 is a diagram illustrating the structure of a 
broadband engine (BE) in accordance with the present 
invention; 

Figure 4 is a diagram illustrating the structure of an 
synergistic processing unit (SPU) in accordance with the 
present invention; 

Figure 5 is a diagram illustrating the structure of a 
processing unit, visualizer (VS) and an optical interface 
in accordance with the present invention; 

Figure 6 is a diagram illustrating one combination of 
processing units in accordance with the present invention; 

Figure 7 illustrates another combination of processing 
units in accordance with the present invention; 



Docket No. AUS920030703US1 6 Atty. Ref. No. roM-3064 

Figure 8 illustrates yet another combination of 
processing units in accordance with the present invention; 

Figure 9 illustrates yet another combination of 
processing units in accordance with the present invention; 

Figure 10 illustrates yet another combination of 
processing units in accordance with the present invention; 

Figure llA illustrates the integration of optical 
interfaces within a chip package in accordance with the 
present invention; 

Figure IIB is a diagram of one configuration of 
processors using the optical interfaces of Figure llA; 

Figure IIC is a diagram of another configuration of 
processors using the optical interfaces of Figure llA; 

Figure 12A illustrates the structure of a memory 
system in accordance with the present invention; 

Figure 12B illustrates the writing of data from a 
first broadband engine to a second broadband engine in 
accordance with the present invention; 

Figure 13 is a diagram of the structure of a shared 
memory for a processing unit in accordance with the present 
invention; 

Figure 14A illustrates one structure for a bank of the 
memory shown in Figur 13; 
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Figure 14B illustrates another structure for a bank of 
the memory shown in Figure 13; 

Figure 15 illustrates a structure for a direct memory 
access controller in accordance with the present invention; 

Figure 16 illustrates an alternative structure for a 
direct memory access controller in accordance with the 
present invention; 

Figures 17-31 illustrate the operation of data 
synchronization in accordance with the present invention; 

Figure 32 is a three-state memory diagram illustrating 
the various states of a memory location in accordance with 
the data synchronization scheme of the-present invention; 

Figure 33 illustrates the structure of a key control 
table for a hardware sandbox in accordance with the present 
invention; 

Figure 34 illustrates a scheme for storing memory 
access keys for a hardware sandbox in accordance with the 
present invention; 

Figure 35 illustrates the structure of a memory access 
control table for a hardware sandbox in accordance with the 
present invention; 

Figure 36 is a flow diagram of the steps for accessing 
a memory sandbox using the key control table of Figure 33 
and the memory access control table of Figure 35; 
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Figure 37 illustrates the structure of a software cell 
in accordance with the present invention; 

Figure 38 is a flow diagr^am of the steps for issuing 
remote procedure calls to SPUs in accordance with the 
present invention; 

Figure 39 illustrates the structure of a dedicated 
pipeline for processing streaming data in accordance with 
the present invention; 

Figure 40 is a flow diagram of the steps performed by 
the dedicated pipeline of Figure 39 in the processing of 
streaming data in accordance with the present invention; 

Figure 41 illustrates an alternative structure for a 
dedicated pipeline for the processing of streaming data in 
accordance with the present invention; 

Figure 42 illustrates a scheme for an absolute timer 
for coordinating the parallel processing of applications 
and data by SPUs in accordance with the present invention; 

Figure 43 is a block diagram illustrating a processing 
element having a main processor and a plurality of 
secondary processors sharing a system memory; 

Figure 44 is a block diagram illustrating a processing 
element having a main processor and a plurality of 
secondary processors sharing a system memory; 

Figure 45 is a flowchart illustrating a method for 
loading data from the disk to the common system memory; 



Docket No. AUS920030703US1 9 Atty. Ref. No. IBM-3064 

Figure 46 is a flowchart illustrating a process for 
parallel processing data in a common system memory with a 
plurality of processors. 

Figure 47 is a block diagram illustrating creation, 
from a system of linear equations, of an equivalent 
augmented matrix; 

Figure 48 is a block diagram illustrating division of 
the linear equations coefficients into blocks and the 
loading of the blocks into a common memory; 

Figure 49 is a table illustrating the matrix 
operations and whether, for a given block, an SPU has 
completed the operation; 

Figure 50 is a block diagram illustrating the SPU's 
accessing of the common memory and performing the matrix 
operations for a given block; 

Figure 51 is a flowchart illustrating the receiving of 
the linear equations coefficients and the loading of the 
coefficients into the common memory; 

Figure 52 is a flowchart illustrating the PU 
determining a set of matrix operations to solve the linear 
differential equations; and 

Figure 53 is a flowchart illustrating the SPUs 
performing the matrix operations on a block-by-block basis. 
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DETAILED DESCRIPTION 

The following is intended to provide a detailed 
description of an example of the invention and should not 
be taken to be limiting of the invention itself. Rather, 
any number of variations may fall within the scope of the 
invention defined in the claims following the description. 

The overall architecture for a computer system 101 in 
accordance with the present invention is shown in Figure 1. 

As illustrated in this figure, system 101 includes 

network 104 to which is connected a plurality of computers 
and computing devices. Network 104 can be a LAN, a global 
network, such as the Internet, or any other computer 
network . 

The computers and computing devices connected to 
network 104 (the network's "members") include, e.g., client 
computers 106, server computers 108, personal digital 
assistants (PDAs) 110, digital television (DTV) 112 and 
other wired or wireless computers and computing devices. 
The processors employed by the members of network 104 are 
constructed from the same common computing module. These 
processors also preferably all have the same ISA and 
perform processing in accordance with the same instruction 
set. The nximber of modules included within any particular 
processor depends upon the processing power required by 
that processor. 

For example, since servers 108 of system 101 perform 
more processing of data and applications than clients 106, 
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servers 108 contain more computing modules than clients 
106. PDAS 110, on the other hand, perform the least amount 
of processing. PDAs 110, therefore, contain the smallest 
number of computing modules. DTV 112 performs a level of 
processing between that of clients 106 and servers 108. DTV 
112, therefore, contains a number of computing modules 
between that of clients 106 and servers 108. As discussed 
below, each computing module contains a processing 
controller and a plurality of identical processing units 
for performing parallel processing of the data and 
applications transmitted over network 104. 

This homogeneous configuration for system 101 
facilitates adaptability, processing speed and processing 
efficiency. Because each member of system 101 performs 
processing using one or more (or some fraction) of the same 
computing module, the particular computer or computing 
device performing the actual processing of data and 
applications is unimportant. The processing of a particular 
application and data, moreover, can be shared among the 
network's members. By uniquely identifying the cells 
comprising the data and applications processed by system 
101 throughout the system, the processing results can be 
transmitted to the computer or computing device requesting 
the processing regardless of where this processing 
occurred. Because the modules performing this processing 
have a common structure and employ a common ISA, the 
computational burdens of an added layer of software to 
achieve compatibility among the processors is avoided. This 
architecture and programming model facilitates the 
processing speed necessary to execute, e.g., real-time, 
multimedia applications. 
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To take further advantage of the processing speeds and 
efficiencies facilitated by system 101, the data and 
applications processed by this system are packaged into 
uniquely identified, uniformly formatted software cells 
102. Each software cell 102 contains, or can contain, both 
applications and data. Each software cell also contains an 
ID to globally identify the cell throughout network 104 and 
system 101. This uniformity of structure for the software 
cells, and the software cells' unique identification 
throughout the network, facilitates the processing of 
applications and data on any computer or computing device 
of the network. For example, a client 106 may formulate a 
software cell 102 but, because of the limited processing 
capabilities of client 106, transmit this software cell to 
a server 108 for processing. Software cells can migrate, 
therefore, throughout network 104 for processing on the 
basis of the availability of processing resources on the 
network . 

The homogeneous structure of processors and software 
cells of system 101 also avoids many of the problems of 
today's heterogeneous networks. For example, inefficient 
programming models which seek to permit processing of 
applications on any ISA using any instruction set, e.g., 
virtual machines such as the Java virtual machine, are 
avoided. System 101, therefore, can implement broadband 
processing far more effectively and efficiently than 
today ' s networks . 

The basic processing module for all members of network 
104 is the processing unit (PU) . Figure 2 illustrates the 
structure of a PU. As shown in this figure, PU 201 
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comprises a processing unit (PU) 203, a direct memory 
access controller (DMAC) 205 and a plurality of synergistic 
processing units (SPUs) , namely, SPU 207, SPU 209, SPU 211, 
SPU 213, SPU 215, SPU 217, SPU 219 and SPU 221* A local PU 
bus 223 transmits data and applications among the SPUs, 
DMAC 205 and PU 203. Local PU bus 223 can have, e.g., a 
conventional architecture or be implemented as a packet 
switch network. Implementation as a packet switch network, 
while requiring more hardware, increases available 
bandwidth. 

PU 201 can be constructed using various methods for 
implementing digital logic. PU 201 preferably is 
constructed, however, as a single integrated circuit 
employing a complementary metal oxide semiconductor (CMOS) 
on a silicon substrate. Alternative materials for 
substrates include gallium arsinide, gallium aluminum 
arsinide and other so-called III-B compounds employing a 
wide variety of dopants. PU 201 also could be implemented 
using superconducting material, e.g., rapid single-flux- 
quantum ( RSFQ ) logic . 

PU 201 is closely associated with a dynamic random 
access memory (DRAM) 225 through a high bandwidth memory 
connection 227. DRAM 225 functions as the main memory for 
PU 201. Although a DRAM 225 preferably is a dynamic random 
access memory, DRAM 225 could be implemented using other 
means, e.g., as a static random access memory (SRAM), a 
magnetic random access memory (MRAM) , an optical memory or 
a holographic memory. DMAC 205 facilitates the transfer of 
data between DRAM 225 and the SPUs and PU of PU 201. As 
further discussed below, DMAC 205 designates for each SPU 
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an exclusive area in DRAM 225 into which only the SPU can 
write data and from which only the SPU can read data. This 
exclusive area is designated a "sandbox." 

PU 203 can be, e.g., a standard processor capable of 
stand-alone processing of data and applications. In 
operation, PU 203 schedules and orchestrates the processing 
of data and applications by the SPUs. The SPUs preferably 
are single instruction, multiple data (SIMD) processors. 
Under the control of PU 2 03, the SPUs perform the 
processing of these data and applications in a parallel and 
independent manner. DMAC 205 controls accesses by PU 203 
and the SPUs to the data and applications stored in the 
shared DRAM 225. Although PU 201 preferably includes eight 
SPUs, a greater or lesser number of SPUs can be employed in 
a PU depending upon the processing power required. Also, a 
number of PUs, such as PU 201, may be joined or packaged 
together to provide enhanced processing power. 

For example, as shown in Figure 3, four PUs may be 
packaged or joined together, e.g., within one or more chip 
packages, to form a single processor for a member of 
network 104. This configuration is designated a broadband 
engine (BE). As shown in Figure 3, BE 301 contains four 
PUs, namely, PU 303, PU 305, PU 307 and PU 309. 
Communications among these PUs are over BE bus 311. Broad 
bandwidth memory connection 313 provides communication 
between shared DRAM 315 and these PUs. In lieu of BE bus 
311, communications among the PUs of BE 301 can occur 
through DRAM' 315 and this memory connection. 
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Input/output (I/O) interface 317 and external bus 319 
provide communications between broadband engine 301 and the 
other members of network 104. Each PU of BE 301 performs 
processing of data and applications in a parallel and 
independent manner analogous to the parallel and 
independent processing of applications and data performed 
by the SPUs of a PU. 

Figure 4 illustrates the structure of an SPU. SPU 402 
includes local memory 406, registers 410, four floating 
point units 412 and four integer units 414. Again, however, 
depending upon the processing power required, a greater or 
lesser number of floating points units 412 and integer 
units 414 can be employed. In a preferred embodiment, local 
memory 406 contains 128 kilobytes of storage, and the 
capacity of registers 410 is 128 . times . 128 bits. Floating 
point units 412 preferably operate at a speed of 32 billion 
floating point operations per second (32 GFLOPS) , and 
integer units 414 preferably operate at a -speed of 32 
billion operations per second (32 GOPS) . 

Local memory 406 is not a cache memory. Local memory 
406 is preferably constructed as an SRAM. Cache coherency 
support for an SPU is unnecessary. A PU may require cache 
coherency support for direct memory accesses initiated by 
the PU. Cache coherency support is not required, however, 
for direct memory accesses initiated by an SPU or for 
accesses from and to external devices. 

SPU 402 further includes bus 404 for transmitting 
applications and data to and from the SPU. In a preferred 
embodiment, this bus is 1,024 bits wide. SPU 402 further 
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includes internal busses 408, 420 and 418. In a preferred 
embodiment, bus 408 has a width of 256 bits and provides 
communications between local memory 406 and registers 410. 
Busses 420 and 418 provide communications between, 
respectively, registers 410 and floating point units 412, 
and registers 410 and integer units 414. In a preferred 
embodiment, the width of busses 418 and 420 from registers 
410 to the floating point or integer units is 384 bits, and 
the width of busses 418 and 420 from the floating point or 
integer units to registers 410 is 128 bits. The larger 
width of these busses from registers 410 to the floating 
point or integer units than from these units to registers 
410 accommodates the larger data flow from registers 410 
during processing. A maximum of three words are needed for 
each calculation. The result of each calculation, however, 
normally is only one word. 

Figures. 5-10 further illustrate the modular structure 
of the processors of the members of network 104. For 
example, as shown in Figure 5, a processor may comprise a 
single PU 502. As discussed above, this PU typically 
comprises a PU, DMAC and eight SPUs. Each SPU includes 
local storage (LS) . On the other hand, a processor may 
comprise the structure of visualizer (VS) 505. As shown in 
Figure 5, VS 505 comprises PU 512, DMAC 514 and four SPUs, 
namely, SPU 516, SPU 518, SPU 520 and SPU 522. The space 
within the chip package normally occupied by the other four 
SPUS of a PU is occupied in this case by pixel engine 508, 
image cache 510 and cathode ray tube controller (CRTC) 504. 
Depending upon the speed of communications required for PU 
502 or VS 505, optical interface 506 also may be included 
on the chip package. 
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using this standardized, modular structure, numerous 
other variations of processors can be constructed easily 
and efficiently. For example, the processor shown in Figure 
6 comprises two chip packages, namely, chip package 602 
comprising a BE and chip package 604 comprising four VSs. 
input/output (I/O) 606 provides an interface between the BE 
of chip package 602 and network 104. Bus 608 provides 
communications between chip package 602 and chip package 
604. input output processor (lOP) 610 controls the flow of 
data into and out of I/O 606. I/O 606 may be fabricated as 
an application specific integrated circuit (ASIC) . The 
output from the VSs is video signal 612. 

Figure 7 illustrates a chip package for a BE 702 with 
two optical interfaces 704 and 706 for providing ultra high 
speed communications to the other members of network 104 
(or other chip packages locally connected) . BE 702 can 
function as, e.g.,- a server on network 104. 

The chip package of Figure 8 comprises two PUs 802 and 
804 and two VSs 806 and 808. An I/O 810 provides an 
interface between the chip package and network 104. The 
output from the chip package is a video signal. This 
configuration may function as, e.g., a graphics work 
station. 

Figure 9 illustrates yet another configuration. This 
configuration contains one-half of the processing power of 
the configuration illustrated in Figure 8 . Instead of two 
PUS, one PU 902 is provided, and instead of two VSs, one VS 
904' is provided. I/O 906 has one-half the bandwidth of the 
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I/O illustrated in Figure 8. Such a processor also may 
function, however, as a graphics work station. 

A final configuration is shown in Figure 10. This 
processor consists of only a single VS 1002 and an I/O 
1004. This configuration may function as, e.g., a PDA. 

Figure llA illustrates the integration of optical 
interfaces into a chip package of a processor of network 
104. These optical interfaces convert optical signals to 
electrical signals and electrical signals to optical 
signals and can be constructed from a variety of materials 
including, e.g., gallium arsinide, aluminum gallium 
arsinide, germanium and other elements or compounds. As 
shown in this figure, optical interfaces 1104 and 1106 are 
fabricated on the chip package of BE 1102. BE bus 1108 
provides communication among the PUs of BE 1102, namely, PU 
1110, PU 1112, PU 1114, PU 1116, and these optical 
interfaces. Optical interface 1104 includes two ports, 
namely, port 1118 and port 1120, and optical interface 1106 
also includes two ports, namely, port 1122 and port 1124. 
Pdrts 1118, 1120, 1122 and 1124 are connected to, 
respectively, optical wave guides 1126, 1128, 1130 and 
1132 . Optical signals are transmitted to and from BE 1102 
through these optical wave guides via the ports of optical 
interfaces 1104 and 1106. 

plurality of BEs can be connected together in various 
configurations using such optical wave guides and the four 
optical ports of each BE. For example, as shown in Figure 
IIB, two or more BEs, e.g., BE 1152, BE 1154 and BE 1156, 
can be connected serially through such optical ports. In 
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this example, optical interface 1166 of BE 1152 is 
connected through its optical ports to the optical ports of 
optical interface 1160 of BE 1154. In a similar manner, the 
optical ports of optical interface 1162 on BE 1154 are 
connected to the optical ports of optical interface 1164 of 
BE 1156. 

A matrix configuration is illustrated in Figure IIC. 
In this configuration, the optical interface of each BE is 
connected to two other EEs. As shown in this figure, one of 
the optical ports of optical interface 1188 of BE 1172 is 
connected to an optical port of optical interface 1182 of 
BE 1176. The other optical port of optical interface 1188 
is connected to an optical port of optical interface 1184 
of BE 1178. In a similar manner, one optical port of 
optical interface 1190 of BE 1174 is connected to the other 
optical port of optical interface 1184 of BE 1178. The 
other optical port of optical interface 1190 is connected 
to an optical port of optical interface 1186 of BE 1180. 
This matrix configuration can be extended in a similar 
manner to other BEs. 

Using either a serial configuration or a matrix 
configuration, a processor for network 104 can be 
constructed of any desired size and power. Of course, 
additional ports can be added to the optical interfaces of 
the BES, or to processors having a greater or lesser number 
of PUS than a BE, to form other configurations. 

Figure 12A illustrates the control system and 
structure for the DRAM of a BE. A similar control system 
and structure is employed in processors having other sizes 



Docket No. AUS920030703US 1 20 



Atty. Ref. No. ffiM-3064 



and containing more or less PUs. As shown in this figure, a 
cross-bar switch connects each DMAC 1210 of the four PUs 
comprising BE 1201 to eight bank controls 1206. Each bank 
control 1206 controls eight banks 1208 (only four are shown 
in the figure) of DRAM 1204. DRAM 1204, therefore, 
comprises a total of sixty- four banks. In a preferred 
embodiment, DRAM 1204 has a capacity of 64 megabytes, and 
each bank has a capacity of 1 megabyte. The smallest 
addressable unit within each bank, in this preferred 
embodiment, is a block of 1024 bits. 

BE 1201 also includes switch unit 1212. Switch unit 
1212 enables other SPUs on BEs closely coupled to BE 1201 
to access DRAM 1204. A second BE, therefore, can be closely 
coupled to a first BE, and each SPU of each BE can address 
twice the niimber of memory locations normally accessible to 
an SPU. The direct reading or writing of data from or to 
the DRAM of a first BE from or to the DRAM of a second BE 
can occur through a switch unit such as switch unit 1212. 

For example, as shown in Figure 12B, to accomplish 
such writing, the SPU of a first BE, e.g., SPU 1220 of BE 
1222, issues a write command to a memory location of a DRAM 
of a second BE, e.g., DRAM 1228 of BE 1226 (rather than, as 
in the usual case, to DRAM 1224 of BE 1222) . DMAC 1230 of 
BE 1222 sends the write command through cross-bar switch 
1221 to bank control 1234, and bank control 1234 transmits 
the command to an external port 1232 connected to bank 
control 1234. DMAC 1238 of BE 1226 receives the write 
command and transfers this command to switch unit 1240 of 
BE 1226. Switch unit 1240 identifies the DRAM address 
contained in the write command and sends the data for 



Docket No. AUS920030703US 1 21 



Atty. Ref. No. ffiM-3064 



storage in this address through bank control 1242 of BE 
1226 to bank 1244 of DRAM 1228. Switch unit 1240, 
therefore, enables both DRAM 1224 and DRAM 1228 to function 
as a single memory space for the SPUs of BE 1226. 

Figure 13 shows the configuration of the sixty-four 
banks of a DRAM. These banks are arranged into eight rows, . 
namely, rows 1302, 1304, 1306, 1308, 1310, 1312, 1314 and 
1316 and eight columns, namely, columns 1320, 1322, 1324, 
1326, 1328, 1330, 1332 and 1334. Each row is controlled by 
a bank controller. Each bank controller, therefore, 
controls eight megabytes of memory. 

Figures. 14A and 14B illustrate different 
configurations for storing and accessing the smallest 
addressable memory unit of a DRAM, e.g., a block of 1024 
bits. In Figure 14A, DMAC 1402 stores in a single bank 1404 
eight 1024 bit blocks 1406. In Figure 14B, on the other 
hand, while DMAC 1412 reads and writes blocks of data 
containing 1024 bits, these blocks are interleaved between 
two banks, namely, bank 1414 and bank 1416. Each of these 
banks, therefore, contains sixteen blocks of data, and each 
block of data contains 512 bits. This interleaving can 
facilitate faster accessing of the DRAM and is useful in 
the processing of certain applications. 

Figure 15 illustrates the architecture for a DMAC 1504 
within a PU. As illustrated in this figure, the structural 
hardware comprising DMAC 1506 is distributed throughout the 
PU such that each SPU 1502 has direct access to a 
structural node 1504 of DMAC 1506. Each node executes the 
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logic appropriate for memory accesses by the SPU to which 
the node has direct access. 

Figure 16 shows an alternative embodiment of the DMAC, 
namely, a non-distributed architecture. In this case, the 
5 structural hardware of DMAC 1606 is centralized. SPUs 1602 

and PU 1604 communicate with DMAC 1606 via local PU bus 
1607. DMAC 1606 is connected through a cross-bar switch to 
a bus 1608. Bus 1608 is connected to DRAM 1610. 

As discussed above, all of the multiple SPUs of a PU 

10 can independently access data in the shared DRAM. As a 

result, a first SPU could be operating upon particular data 
in its local storage at a time during which a second SPU 
requests these data. If the data were provided to the 
second SPU at that time from the shared DRAM, the data 

15 could be invalid because of the first SPU's ongoing 

processing which could change the data's value. If the 
second processor received the data from the shared DRAM at 
that time, therefore, the second processor could generate 
an erroneous result. For example, the data could be a 

20 specific value for a global variable. If the first 

processor changed that value during its processing, the 
second processor would receive an outdated value. A scheme 
is necessary, therefore, to synchronize the SPUs' reading 
and writing of data from and to memory locations within the 

2 5 shared DRAM. This scheme must prevent the reading of data 

from a memory location upon which another SPU currently is 
operating in its local storage and, therefore, which are 
not current, and the writing of data into a memory location 
storing current data. 
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To overcome these problems, for each addressable 
memory location of the DRAM, an additional segment of 
memory is allocated in the DRAM for storing status 
information relating to the data stored in the memory 
5 location. This status information includes a full/empty 

(F/E) bit, the identification of an SPU (SPU ID) requesting 
data from the memory location and the address of the SPU's 
local storage (LS address) to which the requested data 
should be read. An addressable memory location of the DRAM 
10 can be of any size. In a preferred embodiment, this size is 

1024 bits. 

The setting of the F/E bit to 1 indicates that the 
data stored in the associated memory location are current. 
The setting of the F/E bit to 0, on the other hand, 

15 indicates that the data stored in the associated memory 

location are not current. If an SPU requests the data when 
this bit is set to 0, the SPU is prevented from immediately 
reading the data. In this case, an SPU ID identifying the 
SPU requesting the data, and an LS address identifying the 

20 memory location within the local storage of this SPU to 

which the data are to be read when the data become current, 
are entered int.o the additional memory segment. 

An additional memory segment also is allocated for 
each memory location within the local storage of the SPUs. 

25 This additional memory segment stores one bit, designated 

the "busy bit." The busy bit is used to reserve the 
associated LS memory location for the storage of specific 
data to be retrieved from the DRAM. If the busy bit is set 
to 1 for a particular memory location in local storage, the 

30 SPU can use this memory location only for the writing of 
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these specific data. On the other hand, if the busy bit is 
set to 0 for a particular memory location in local storage, 
the SPU can use this memory location for the writing of any 
data. 

Examples of the manner in which the F/E bit, the SPU 
ID, the LS address and the busy bit are used to synchronize 
the reading and writing of data from and to the shared DRAM 
of a PU are illustrated in Figures. 17-31. 

As shown in Figure 17, one or more PUs, e.g., PU 1720, 
interact with DRAM 1702. PU 1720 includes SPU 1722 and SPU 
1740. SPU 1722 includes control logic 1724, and SPU 1740 
includes control logic 1742. SPU 1722 also includes local 
storage 1726. This local storage includes a plurality of 
addressable memory locations 1728. SPU 1740 includes local 
storage 1744, and this local storage also includes a 
plurality of addressable memory locations 1746. All of 
these addressable memory locations preferably are 1024 bits 
in size. 

An additional segment of memory is associated with 
each LS addressable memory location. For example, memory 
segments 1729 and 1734 are associated with, respectively, 
local memory locations 1731 and 1732, and memory segment 
1752 is associated with local memory location 1750. A "busy 
bit," as discussed above, is stored in each of these 
additional memory segments. Local memory location 1732 is 
shown with several Xs to indicate that this location 
contains data. 
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DRAM 1702 contains a plurality of addressable memory 
locations 1704, including memory locations 1706 and 1708. 
These memory locations preferably also are 1024 bits in 
size. An additional segment of memory also is associated 
with each of these memory locations. For example, 
additional memory segment 1760 is associated with memory 
location 1706, and additional memory segment 1762 is 
associated with memory location 1708. Status information 
relating to the data stored in each memory location is 
stored in the memory segment associated with the memory 
location. This status information includes, as discussed 
above, the F/E bit, the SPU ID and the LS address. For 
example, for memory location 1708, this status information 
includes F/E bit 1712, SPU ID 1714 and LS address 1716. 

Using the status information and the busy bit, the 
synchronized reading and writing of data from and to the 
shared DRAM among the SPUs of a PU, or a group of PUs, can 
be achieved. 

Figure 18 illustrates the initiation of the 
synchronized writing of data from LS memory location 1732 
of SPU 1722 to memory location 1708 of DRAM 1702. Control 
1724 of SPU 1722 initiates the synchronized writing of 
these data. Since memory location 1708 is empty, F/E bit 
1712 is set to 0. As a result, the data in LS location 1732 
can be written into memory location 1708. If this bit were 
set to 1 to indicate that memory location 1708 is full and 
contains current, valid data, on the other hand, control 
1722 would receive an error message and be prohibited from 
writing data into this memory location. 
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The result of the successful synchronized writing of 
the data into memory location 1708 is shown in Figure 19. 
The written data are stored in memory location 1708, and 
F/E bit 1712 is set to 1. This setting indicates that 
memory location 1708 is full and that the data in this 
memory location are current and valid. 

Figure 20 illustrates the initiation of the 
synchronized reading of data from memory location 1708 of 
DRAM 1702 to LS memory location 1750 of local storage 1744. 
To initiate this reading, the busy bit in memory segment 
1752 of LS memory location 1750 is set to 1 to reserve this 
memory location for these data. The setting of this busy 
bit to 1 prevents SPU 1740 from storing other data in this 
memory location. 

As shown in Figure 21, control logic 1742 next issues 
a synchronize read command for memory location 1708 of DRAM 
1702. Since F/E bit 1712 associated with this memory 
location is set to 1, the data stored in memory location 
1708 are considered current and valid. As a result, in 
preparation for transferring the data from memory location 
1708 to LS memory location 1750, F/E bit 1712 is set to 0. 
This setting is shown in Figure 22. The setting of this bit 
to 0 indicates that, following the reading of these data, 
the data in memory location 1708 will be invalid. 

As shown in Figure 23, the data within memory location 
1708 next are read from memory location 1708 to LS memory 
location 1750. Figure 24 shows the final state. A copy of 
the data in memory location 1708 is stored in LS memory 
location 1750. F/E bit 1712 is set to 0 to indicate that 
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the data in memory location 1708 are invalid. This 
invalidity is the result of alterations to these data to be 
made by SPU 1740. The busy bit in memory segment 1752 also 
is set to 0. This setting indicates that LS memory location 
1750 now is available to SPU 1740 for any purpose, i.e., 
this LS memory location no longer is in a reserved state 
waiting for the receipt of specific data. LS memory 
location 1750, therefore, now can be accessed by SPU 1740 
for any purpose. 

Figures. 25-31 illustrate the synchronized reading of 
data from a memory location of DRAM 1702, e.g., memory 
location 1708, to an LS memory location of an SPU's local 
storage, e.g., LS memory location 1752 of local storage 
1744, when the F/E bit for the memory location of DRAM 1702 
is set to 0 to indicate that the data in this memory 
location are not current or valid. As shown in Figure 25, 
to initiate this transfer, the busy bit in memory segment 
1752 of LS memory location 1750 is set to 1 to reserve this 
LS memory location for this transfer of data. As shown in 
Figure 26, control logic 1742 next issues a synchronize 
read command for memory locatiori 1708 of DRAM 1702. Since 
the F/E bit associated with this memory location, F/E bit 
1712, is set to 0, the data stored in memory location 1708 
are invalid. As a result, a signal is transmitted to 
control logic 1742 to block the immediate reading of data 
from this memory location. 

AS shown in Figure 27, the SPU ID 1714 and LS address 
1716 for this read command next are written into memory 
segment 1762. In this case, the SPU ID for SPU 1740 and the 
LS memory location for LS memory location 1750 are written 
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into memory segment 1762. When the data within memory 
location 1708 become current, therefore, this SPU ID and LS 
memory location are used for determining the location to 
which the current data are to be transmitted. 

The data in memory location 1708 become valid and 
current when an SPU writes data into this memory location. 
The synchronized writing of data into memory location 1708 
from, e.g., memory location 1732 of SPU 1722, is 
illustrated in Figure 28. This synchronized writing of 
these data is permitted because F/E bit 1712 for this 
memory location is set to 0. 

AS shown in Figure 29, following this writing, the 
data in memory location 1708 become current and valid. SPU 
ID 1714 and LS address 1716 from memory segment 1762, 
therefore, immediately are read from memory segment 1762, 
and this information then is deleted from this segment. F/E 
bit 1712 also is set to 0 in anticipation of the immediate 
reading of the data in memory location 1708. As shown in 
Figure 30, upon reading SPU ID 1714 and LS address 1716, 
this information immediately is used for reading the valid 
data in memory location 1708 to LS memory location 1750 of 
SPU 1740. The final state is shown in Figure 31. This 
figure shows the valid data from memory location 1708 
copied to memory location 1750, the busy bit in memory 
segment 1752 set to 0 and F/E bit 1712 in memory segment 
1762 set to 0. The setting of this busy bit to 0 enables LS 
memory location 1750 now to be accessed by SPU 1740 for any 
purpose. The setting of this F/E bit to 0 indicates that 
the data in memory location 1708 no longer are current and 
valid. 
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Figure 32 summarizes the operations described above 
and the various states of a memory location of the DRAM 
based upon the states of the F/E bit, the SPU ID and the LS 
address stored in the memory segment corresponding to the 
memory location. The memory location can have three states. 
These three states are an empty state 3280 in which the F/E 
bit is set to 0 and no information is provided for the SPU 
ID or the LS address, a full state 3282 in which the F/E 
bit is set to 1 and no information is provided for the SPU 
ID or LS address and a blocking state 3284 in which the F/E 
bit is set to 0 and information is provided for the SPU ID 
and LS address. 

As shown in this figure, in empty state 3280, a 
synchronized writing operation is permitted and results in 
a transition to full state 3282. A synchronized reading 
operation, however, results in a transition to the blocking 
state 3284 because the data in the memory location, when 
the memory location is in the empty state, are not current. 

In full state 3282, a synchronized reading operation 
is permitted and results in a transition to empty state 
3280. On the other hand, a synchronized writing operation 
in full state 3282 is prohibited to prevent overwriting of 
valid data. If such a writing operation is attempted in 
this state, no state change occurs and an error message is 
transmitted to the SPU's corresponding control logic. 

In blocking state 3284, the synchronized writing of 
data into the memory location is permitted and results in a 
transition to empty state 3280. On the other hand, a 
synchronized reading operation in blocking state 3284 is 
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prohibited to prevent a conflict with the earlier 
synchronized reading operation which resulted in this 
state. If a synchronized reading operation is attempted in 
blocking state 3284, no state change occurs and an error 
5 message is transmitted to the SPU's corresponding control 
logic . 

The scheme described above for the synchronized 
reading and writing of data from and to the shared DRAM 
also can be used for eliminating the computational 

10 resources normally dedicated by a processor for reading 

data from, and writing data to, external devices. This 
input/output (I/O) function could be performed by a PU. 
However, using a modification of this synchronization 
scheme, an SPU running an appropriate program can perform 

15 this function. For example, using this scheme, a PU 

receiving an interrupt request for the transmission of data 
from an I/O interface initiated by an external device can 
delegate the handling of this request to this SPU. The SPU 
then issues a synchronize write command to the I/O 

20 interface. This interface in turn signals the external 

device that data now can be written into the DRAM. The SPU 
next issues a synchronize read command to the DRAM to set 
the DRAM'S relevant memory space into a blocking state. The 
SPU also sets to 1 the busy bits for the memory locations 

25 of the SPU's local storage needed to receive the data. In 

the blocking state, the additional memory segments 
associated with the DRAM's relevant memory space contain 
the SPU's ID and the address of the relevant memory 
locations of the SPU's local storage. The external device 

30 next issues a synchronize write command to write the data 

directly to the DRAM's relevant memory space. Since this 
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memory space is in the blocking state, the data are 
immediately read out of this space into the memory 
locations of the SPU's local storage identified in the 
additional memory segments. The busy bits for these memory 
locations then are set to 0. When the external device 
completes writing of the data, the SPU issues a signal to 
the PU that the transmission is complete. 

using this scheme, therefore, data transfers from 
external devices can be processed with minimal 
computational load on the PU. The SPU delegated this 
function, however, should be able to issue an interrupt 
request to the PU, and the external device should have 
direct access to the DRAM. 

The DRAM of each PU includes a plurality of 
"sandboxes." A sandbox defines an area of the shared DRAM 
beyond which a particular SPU, or set of SPUs, cannot read 
or write data. These sandboxes provide security against the 
corruption of data being processed by one SPU by data being 
processed by another SPU. These sandboxes also permit the 
downloading of software cells from network 104 into a 
particular sandbox without the possibility of the software 
cell corrupting data throughout the DRAM. In the present 
invention, the sandboxes are implemented in the hardware of 
the DRAMS and DMACs. By implementing these sandboxes in 
this hardware rather than in software, advantages in speed 
and security are obtained. 

The PU of a PU controls the sandboxes assigned to the 
SPUS. Since the PU normally operates only trusted programs, 
such as an operating system, this scheme does not 
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jeopardize security. In accordance with this scheme, the PU 
builds and maintains a key control table. This key control 
table is illustrated in Figure 33. As shown in this figure, 
each entry in key control table 3302 contains an 
identification (ID) 3304 for an SPU, an SPU key 3306 for 
that SPU and a key mask 3308. The use of this key mask is 
explained below. Key control table 3302 preferably is 
stored in a relatively fast memory, such as a static random 
access memory (SRAM), and is associated with the DMAC. The 
entries in key control table 3302 are controlled by the PU. 
When an SPU requests the writing of data to, or the reading 
of data from, a particular storage location of the DRAM, 
the DMAC evaluates the SPU key 3306 assigned to that SPU in 
key control table 3302 against a memory access key 
associated with that storage location. 

As shown in Figure 34, a dedicated memory segment 3410 
is assigned to each addressable storage location 3406 of a 
DRAM 3402. A memory access key 3412 for the storage 
location is stored in this dedicated memory segment. As 
discussed above, a further additional dedicated memory 
segment 3408, also associated with each addressable storage 
location 3406, stores synchronization information for 
writing data to, and reading data from, the storage- 
location. 

In operation, an SPU issues a DMA command to the DMAC. 
This command includes the address of a storage location 
3406 of DRAM 3402. Before executing this command, the DMAC 
looks up the requesting SPU's key 3306 in key control table 
3302 using the SPU's ID 3304. The DMAC then compares the 
SPU key 3306 of the requesting SPU to the memory access key 
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3412 stored in the dedicated memory segment 3410 associated 
with the storage location of the DRAM to which the SPU 
seeks access. If the two keys do not match, the DMA command 
is not executed. On the other hand, if the two keys match, 
the DMA command proceeds and the requested memory access is 
executed . 

An alternative embodiment is illustrated in Figure 35. 
In this embodiment, the PU also maintains a memory access 
control table 3502. Memory access control table 3502 
contains an entry for each sandbox within the DRAM. In the 
particular example of Figure 35, the DRAM contains 64 
sandboxes. Each entry in memory access control table 3502 
contains an identification (ID) 3504 for a sandbox, a base 
memory address 3506, a sandbox size 3508, a memory access 
key 3510 and an access key mask 3512. Base memory address 
3506 provides the address in the DRAM which starts a 
particular memory sandbox. Sandbox size 3508 provides the 
size of the sandbox and, therefore, the endpoint of the 
particular sandbox. 

Figure 36 is a flow diagram of the steps for executing 
a DMA command using key control table 3302 and memory 
access control table 3502.. In step 3602, an SPU issues a 
DMA command to the DMAC for access to a particular memory 
location or locations within a sandbox. This command 
includes a sandbox ID 3504 identifying the particular 
sandbox for which access is requested. In step 3604, the 
DMAC looks up the requesting SPU's key 3306 in key control 
table 3302 using the SPU's ID 3304. In step 3606, the DMAC 
uses the sandbox ID 3504 in the command to look up in 
memory access control table 3502 the memory access key 3510 
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associated with that sandbox. In step 3608, the DMAC 
compares the SPU key 3306 assigned to the requesting SPU to 
the access key 3510 associated with the sandbox. In step 
3610, a determination is made of whether the two keys 
match. If the two keys do not match, the process moves to 
step 3612 where the DMA command does not proceed and an 
error message is sent to either the requesting SPU, the PU 
or both, on the other hand, if at step 3610 the two keys 
are found to match, the process proceeds to step 3614 where 
the DMAC executes the DMA command. 

The key masks for the SPU keys and the memory access 
keys provide greater flexibility to this system. A key mask 
for a key converts a masked bit into a wildcard. For 
example, if the key mask 3308 associated with an SPU key 
3306 has its last two bits set to "mask, " designated by, 
e.g., setting these bits in key mask 3308 to 1, the SPU key 
can be either a 1 or a 0 and still match the memory access 
key. For example, the SPU key might be 1010. This SPU key 
normally allows access only to a sandbox having an access 
key of 1010. If the SPU key mask for this SPU key is set to 
0001, however, then this SPU key can be used to gain access 
to sandboxes having an access key of either 1010 or 1011. 
Similarly, an access key 1010 with a mask set to 0001 can 
be accessed by an SPU with an SPU key of either 1010 or 
1011. Since both the SPU key mask and the memory key mask 
can be used simultaneously, numerous variations of 
accessibility by the SPUs to the sandboxes can be 
established. 



The present invention also provides a new programming 
model for the processors of system 101. This programming 
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model employs software cells 102. These cells can be 
transmitted to any processor on network 104 for processing. 
This new programming model also utilizes the unique modular 
architecture of system 101 and the processors of system 
101. 

Software cells are processed directly by the SPUs from 
the SPU's local storage. The SPUs do not directly operate 
on any data or programs in the DRAM. Data and programs in 
the DRAM are read into the SPU's local storage before the 
SPU processes these data and programs. The SPU's local 
storage, therefore, includes a program counter, stack and 
other software elements for executing these programs. The 
PU controls the SPUs by issuing direct memory access (DMA) 
commands to the DMAC. 

The structure of software cells 102 is illustrated in 
Figure 37. As shown in this figure, a software cell, e.g., 
software cell 3702, contains routing information section 
3704 and body 3706. The information contained in routing 
information section 3704 is dependent upon the protocol of 
network 104. Routing information section 3704 contains 
header 3708, destination ID 3710, source ID 3712 and reply 
ID 3714. The destination ID includes a network address. 
Under the TCP/IP protocol, e.g., the network address is an 
Internet protocol (IP) address. Destination ID 3710 further 
includes the identity of the PU and SPU to which the cell 
should be transmitted for processing. Source ID 3712 
contains a network address and identifies the PU and SPU 
from which the cell originated to enable the destination PU 
and SPU to obtain additional information regarding the cell 
if necessary. Reply ID 3714 contains a network address and 
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identifies the PU and SPU to which queries regarding the 
cell, and the result of processing of the cell, should be 
directed. 

Cell body 3706 contains information independent of the 
network's protocol. The exploded portion of Figure 37 shows 
the details of cell body 3706. Header 3720 of cell body 
3706 identifies the start of the cell body. Cell interface 
3722 contains information necessary for the cell's 
utilization. This information includes global unique ID 
3724, required SPUs 3726, sandbox size 3728 and previous 
cell ID 3730. 

Global unique ID 3724 uniquely identifies software 
cell 3702 throughout network 104. Global unique ID 3724 is 
generated on the basis of source ID 3712, e.g. the unique 
identification of a PU or SPU within source ID 3712, and 
the time and date of generation or transmission of software 
cell 3702. Required SPUs 3726 provides the minimxim number 
of SPUs required to execute the cell. Sandbox size 3728 
provides the amount of protected memory in the required 
SPUs' associated DRAM necessary to execute the cell. 
Previous cell ID 3730 provides the identity of a previous 
cell in a group of cells requiring sequential execution, 
e.g., streaming data. 

implementation section 3732 contains the cell's core 
information. This information includes DMA command list 
3734, programs 3736 and data 3738. Programs 3736 contain 
the programs to be run by the SPUs (called "spulets"), 
e.g., SPU programs 3760 and 3762, and data 3738 contain the 
data to be processed with these programs. DMA command list 
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3734 contains a series of DMA commands needed to start the 
programs. These DMA commands include DMA commands 3740, 
3750, 3755 and 3758. The PU issues these DMA commands to 
the DMAC. 

DMA command 3740 includes VID 3742. VID 3742 is the 
virtual ID of an SPU which is mapped to a physical ID when 
the DMA commands are issued. DMA command 3740 also includes 
load command 3744 and address 3746. Load command 3744 
directs the SPU to read particular information from the 
DRAM into local storage. Address 3746 provides the virtual 
address in the DRAM containing this information. The 
information can be, e.g., programs from programs section 
3736, data from data section 3738 or other data. Finally, 
DMA command 3740 includes local storage address 3748. This 
address identifies the address in local storage where the 
information should be loaded. DMA commands 3750 contain 
similar information. Other DMA commands are also possible. 

DMA command list 3734 also includes a series of kick 
commands, e.g., kick commands 3755 and 3758. Kick commands 
are commands issued by a PU to an SPU to initiate the 
processing of a cell. DMA kick command 3755 includes 
virtual SPU ID 3752, kick command 3754 and program counter 
3756. Virtual SPU ID 3752 identifies the SPU to be kicked, 
kick command 3754 provides the relevant kick command and 
program counter 3756 provides the address for the program 
counter for executing the program. DMA kick command 3758 
provides similar information for the same SPU or another 
SPU. 
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As noted, the PUs treat the SPUs as independent 
processors, not co-processors. To control processing by the 
SPUs, therefore, the PU uses commands analogous to remote 
procedure calls. These commands are designated "SPU Remote 
Procedure Calls" (SRPCs) . A PU implements an SRPC by 
issuing a series of DMA commands to the DMAC. The DMAC 
loads the SPU program and its associated stack frame into 
the local storage of an SPU. The PU then issues an initial 
kick to the SPU to execute the SPU Program. 

Figure 38 illustrates the steps of an SRPC for 
executing an spulet. The steps performed by the PU in 
initiating processing of the spulet by a designated SPU are 
shown in the first portion 3802 of Figure 38, and the steps 
performed by the designated SPU in processing the spulet 
are shown in the second portion 3804 of Figure 38. 

In step 3810, the PU evaluates the spulet and then 
designates an SPU for processing the spulet. In step 3812, 
the PU allocates space in the DRAM for executing the spulet 
by issuing a DMA command to the DMAC to set memory access 
keys for the necessary sandbox or sandboxes. In step 3814, 
the PU enables an interrupt request for the designated SPU 
to signal completion of the spulet. In step 3818, the PU 
issues a DMA command to the DMAC to load the spulet from 
the DRAM to the local storage of the SPU. In step 3820, the 
DMA command is executed, and the spulet is read from the 
DRAM to the SPU's local storage. In step 3822, the PU 
issues a DMA command to the DMAC to load the stack frame 
associated with the spulet from the DRAM to the SPU's local 
storage. In step 3823, the DMA command is executed, and the 
stack frame is read from the DRAM to the SPU's local 
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storage. In step 3824, the PU issues a DMA command for the 
DMAC to assign a key to the SPU to allow the SPU to read 
and write data from and to the hardware sandbox or 
sandboxes designated in step 3812. In step 3826, the DMAC 
updates the key control table (KTAB) with the key assigned 
to the SPU. In step 3828, the PU issues a DMA command 
"kick" to the SPU to start processing of the program. Other 
DMA commands may be issued by the PU in the execution of a 
particular SRPC depending upon the particular spulet. 

As indicated above, second portion 3804 of Figure 38 
illustrates the steps performed by the SPU in executing the 
spulet. In step 3830, the SPU begins to execute the spulet 
in response to the kick command issued at step 3828. In 
step 3832, the SPU, at the direction of the spulet, 
evaluates the spulet ' s associated stack frame. In step 
3834, the SPU issues multiple DMA commands to the DMAC to 
load data designated as needed by the stack frame from the 
DRAM to the SPU's local storage. In step 3836, these DMA 
commands are executed, and the data are read from the DRAM 
to the SPU's local storage. In step 3838, the SPU executes 
the spulet and generates a result. In step 3840, the SPU 
issues a DMA command to the DMAC to store the result in the 
DRAM. In step 3842, the DMA command is executed and the 
result of the spulet is written from the SPU's local 
storage to the DRAM. In step 3844, the SPU issues an 
interrupt request to the PU to signal that the SRPC has 
been completed. 

The ability of SPUs to perform tasks independently 
under the direction of a PU enables a PU to dedicate a 
group of SPUs, and the memory resources associated with a 
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group of SPUs, to performing extended tasks. For example, a 
PU can dedicate one or more SPUs, and a group of memory 
sandboxes associated with these one or more SPUs, to 
receiving data transmitted over network 104 over an 
extended period and to directing the data received during 
this period to one or more other SPUs and their associated 
memory sandboxes for further processing. This ability is 
particularly advantageous to processing streaming data 
transmitted over network 104, e.g., streaming MPEG or 
streaming ATRAC audio or video data. A PU can dedicate one 
or more SPUs and their associated memory sandboxes to 
receiving these data and one or more other SPUs and their 
associated memory sandboxes to decompressing and further 
processing these data. In other words, the PU can establish 
a dedicated pipeline relationship among a group of SPUs and 
their associated memory sandboxes for processing such data. 

In order for such processing to be performed 
efficiently, however, the pipeline's dedicated SPUs and 
memory sandboxes should remain dedicated to the pipeline 
during periods in which processing of spulets comprising 
the data stream does not occur. In other words, the 
dedicated SPUs and their associated sandboxes should be 
placed in a reserved state during these periods. The 
reservation of an SPU and its associated memory sandbox or 
sandboxes upon completion of processing of an spulet is 
called a "resident termination." A resident termination 
occurs in response to an instruction from a PU. 

Figures. 39, 40A and 40B illustrate the establishment 
of a dedicated pipeline structure comprising a group of 
SPUs and their associated sandboxes for the processing of 
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streaming data, e.g., streaming MPEG data. As shown in 
Figure 39, the components of this pipeline structure 
include PU 3902 and DRAM 3918. PU 3902 includes PU 3904, 
DMAC 3906 and a plurality of SPUs, including SPU 3908, SPU 
3910 and SPU 3912. Communications among PU 3904, DMAC 3906 
and these SPUs occur through PU bus 3914. Wide bandwidth 
bus 3916 connects DMAC 3906 to DRAM 3918. DRAM 3918 
includes a plurality of sandboxes, e.g., sandbox 3920, 
sandbox 3922, sandbox 3924 and sandbox 3926. 

Figure 40A illustrates the steps for establishing the 
dedicated pipeline. In step 4010, PU 3904 assigns SPU 3908 
to process a network spulet. A network spulet comprises a 
program for processing the network protocol of network 104. 
In this case, this protocol is the Transmission Control 
Protocol /Internet Protocol (TCP/IP) . TCP/IP data packets 
conforming to this protocol are transmitted over network 
104. Upon receipt, SPU 3908 processes these packets and 
assembles the data in the packets into software cells 102. 
In step 4012, PU 3904 instructs SPU 3908 to perform 
resident terminations upon the completion of the processing 
of the network spulet. In step 4014, PU 3904 assigns PUs 
3910 and 3912 to process MPEG spulets . In step 4015, PU 
3904 instructs SPUs 3910 and 3912 also to perform resident 
terminations upon the completion of the processing of the 
MPEG spulets. In step 4016, PU 3904 designates sandbox 3920 
as a source sandbox for access by SPU 3908 and SPU 3910. In 
step 4018, PU 3904 designates sandbox 3922 as a destination 
sandbox for access by SPU 3910. In step 4020, PU 3904 
designates sandbox 3924 as a source sandbox for access by 
SPU 3908 and SPU 3912. In step 4022, PU 3904 designates 
sandbox 3926 as a destination sandbox for access by SPU 
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3912. In step 4024, SPU 3910 and SPU 3912 send synchronize 
read commands to blocks of memory within, respectively, 
source sandbox 3920 and source sandbox 3924 to set these 
blocks of memory into the blocking state. The process 
finally moves to step 4028 where establishment of the 
dedicated pipeline is complete and the resources dedicated 
to the pipeline are reserved. SPUs 3908, 3910 and 3912 and 
their associated sandboxes 3920, 3922, 3924 and 3926, 
therefore, enter the reserved state. 

Figure 40B illustrates the steps for processing 
streaming MPEG data by this dedicated pipeline. In step 
4030, SPU 3908, which processes the network spulet, 
receives in its local storage TCP/IP data packets from 
network 104. In step 4032, SPU 3908 processes these TCP/IP 
data packets and assembles the data within these packets 
into software cells 102. In step 4034, SPU 3908 examines 
header 3720 (Figure 37) of the software cells to determine 
whether the cells contain MPEG data. If a cell does not 
contain MPEG data, then, in step 4036, SPU 3908 transmits 
the cell to a general purpose sandbox designated within 
DRAM 3918 for processing other data by other SPUs not 
included within the dedicated pipeline. SPU 3908 also 
notifies PU 3904 of this transmission. 

On the other hand, if a software cell contains MPEG 
data, then, in step 4038, SPU 3908 examines previous cell 
ID 3730 (Figure 37) of the cell to identify the MPEG data 
stream to which the cell belongs. In step 4040, SPU 3908 
chooses an SPU of the dedicated pipeline for processing of 
the cell. In this case, SPU 3908 chooses SPU 3910 to 
process these data. This choice is based upon previous cell 
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ID 3730 and load balancing factors. For example, if 
previous cell ID 3730 indicates that the previous software 
cell of the MPEG data stream to which the software cell 
belongs was sent to SPU 3910 for processing, then the 
5 present software cell normally also will be sent to SPU 
3910 for processing. In step 4042, SPU 3908 issues a 
synchronize write command to write the MPEG data to sandbox 
3920. Since this sandbox previously was set to the blocking 
state, the MPEG data, in step 4044, automatically is read 

10 from sandbox 3920 to the local storage of SPU 3910. In step 

4046, SPU 3910 processes the MPEG data in its local storage 
to generate video data. In step 4048, SPU 3910 writes the 
video data to sandbox 3922. In step 4050, SPU 3910 issues a 
synchronize read command to sandbox 3920 to prepare this 

15 sandbox to receive additional MPEG data. In step 4052, SPU 

3910 processes a resident termination. This processing 
causes this SPU to enter the reserved state during which 
the SPU waits to process additional MPEG data in the MPEG 
data stream. 

20 Other dedicated structures can be established among a 

group of SPUs and their associated sandboxes for processing 
other types of data. For example, as shown in Figure 41, a 
dedicated group of SPUs, e.g., SPUs 4102, 4108 and 4114, 
can be established for performing geometric transformations 

25 upon three dimensional objects to generate two dimensional 

display lists. These two dimensional display lists can be 
further processed (rendered) by other SPUs to generate 
pixel data. To perform this processing, sandboxes are 
dedicated to SPUs 4102, 4108 and 4114 for storing the three 

30 dimensional objects and the display lists resulting from 

the processing of these objects. For example, source 
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sandboxes 4104, 4110 and 4116 are dedicated to storing the 
three dimensional objects processed by, respectively, SPU 
4102, SPU 4108 and SPU 4114. In a similar manner, 
destination sandboxes 4106, 4112 and 4118 are dedicated to 
storing the display lists resulting from the processing of 
these three dimensional objects by, respectively, SPU 4102, 
SPU 4108 and SPU 4114. 

Coordinating SPU 4120 is dedicated to receiving in its 
local storage the display lists from destination sandboxes 
4106, 4112 and 4118. SPU 4120 arbitrates among these 
display lists and sends them to other SPUs for the 
rendering of pixel data. 

The processors of system 101 also employ an absolute 
timer. The absolute timer provides a clock signal to the 
SPUs and other elements of a PU which is both independent 
of, and faster than, the clock signal driving these 
elements. The use of this absolute timer is illustrated in 
Figure 42. 

As shown in this figure, the absolute timer 
establishes a time budget for the performance of tasks by 
the SPUs. This time budget provides a time for completing 
these tasks which is longer than that necessary for the 
SPUs' processing of the tasks. As a result, for each task, 
there is, within the time budget, a busy period and a 
standby period. All spulets are written for processing on 
the basis of this time budget regardless of the SPUs ' 
actual processing time or speed. 
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For example, for a particular SPU of a PU, a 
particular task may be performed during busy period 4202 of 
time budget 4204. Since busy period 4202 is less than time 
budget 4204, a standby period 4206 occurs during the time 
budget. During this standby period, the SPU goes into a 
sleep mode during which less power is consumed by the SPU. 

The results of processing a task are not expected by 
other SPUs, or other elements of a PU, until a time budget 
4204 expires. Using the time budget established by the 
absolute timer, therefore, the results of the SPUs' 
processing always are coordinated regardless of the SPUs' 
actual processing speeds. 

In the future, the speed of processing by the SPUs 
will become faster. The time budget established by the 
absolute timer, however, will remain the same. For example, 
as shown in Figure 42, an SPU in the future will execute a 
task in a shorter period and, therefore, will have a longer 
standby period. Busy period 4208, therefore, is shorter 
than busy period 4202, and standby period 4210 is longer 
than standby period 4206. However, since programs are 
written for processing on the basis of the same time budget 
established by the absolute timer, coordination of the 
results of processing among the SPUs is maintained. As a 
result, faster SPUs can process programs written for slower 
SPUs without causing conflicts in the times at which the 
results of this processing are expected. 

In lieu of an absolute timer to establish coordination 
among the SPUs, the PU, or one or more designated SPUs, can 
analyze the particular instructions or microcode being 
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executed by an SPU in processing an spulet for problems in 
the coordination of the SPUs ' parallel processing created 
by enhanced or different operating speeds. "No operation" 
("NOOP") instructions can be inserted into the instructions 
and executed by some of the SPUs to maintain the proper 
sequential completion of processing by the SPUs expected by 
the spulet. By inserting these NOOPs into the instructions, 
the correct timing for the SPUs' execution of all 
instructions can be maintained. 

Figure 43 is a block diagram illustrating a processing 
element having a main processor and a plurality of 
secondary processors sharing a system memory. Processor 
Element (PE) 4305 includes processing unit (PU) 4310, 
which, in one embodiment, acts as the main processor and 
runs an operating system. Processing unit 4310 may be, for 
example, a Power PC core executing a Linux operating 
system, PE 4305 also includes a plurality of synergistic 
processing complex's (SPCs) such as SPCs 4345, 4365, and 
4385. The SPCs include synergistic processing units (SPUs) 
that act as secondary processing units to PU 4310, a memory 
storage unit, and local storage. For example, SPC 4345 
includes SPU 4360, MMU 4355, and local storage 4359; SPC 
4365 includes SPU 4370, MMU 4375, and local storage 4379; 
and SPC 4385 includes SPU 4390, MMU 4395, and local storage 
4399. 

Each SPC may be configured to perform a different 
task, and accordingly, in one embodiment, each SPC may be 
accessed using different instruction sets. If PE 4305 is 
being used in a wireless communications system, for 
example, each SPC may be responsible for separate 
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processing tasks, such as modulation, chip rate processing, 
encoding, network interfacing, etc. In another embodiment, 
the SPCs may have identical instruction sets and may be 
used in parallel with each other to perform operations 
benefiting from parallel processing. 

PE 4305 may also include level 2 cache, such as L2 
cache 4315, for the use of PU 4310. In addition, PE 4305 
includes system memory 4320, which is shared between PU 
4310 and the SPUs. System memory 4320 may store, for 
example, an image of the running operating system (which 
may include the kernel) , device drivers, I/O configuration, 
etc., executing applications, as well as other data. 
System memory 4320 includes the local storage units of one 
or more of the SPCs, which are mapped to a region of system 
memory 4320. For example, local storage 4359 may be mapped 
to mapped region 4335, local storage 4379 may be mapped to 
mapped region 4340, and local storage 4399 may be mapped to 
mapped region 4342. PU 4310 and the SPCs communicate with 
each other and system memory 4320 through bus 4317 that is 
configured to pass data between these devices. 

The MMUs are responsible for transferring data between 
an SPU's local store and the system memory. In one 
embodiment, an MMU includes a direct memory access (DMA) 
controller configured to perform this function. PU 4310 
may program the MMUs to control which memory regions are 
available to each of the MMUs. By changing the mapping 
available to each of the MMUs, the PU may control which SPU 
has access to which region of system memory 4320. In this 
manner, the PU may, for example, designate regions of the 
system memory as private for the exclusive use of a 
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particular SPU. In one embodiment, the SPUs' local stores 
may be accessed by PU 4310 as well as by the other SPUs 
using the memory map. In one embodiment, PU 4310 manages 
the memory map for the common system memory 4320 for all 
the SPUs. The memory map table may include PU 4310 's L2 
Cache 4315, system memory 4320, as well as the SPUs' shared 
local stores. 

In one embodiment, the SPUs process data under the 
control of PU 4310. The SPUs may be, for example, digital 
signal processing cores, microprocessor cores, micro 
controller cores, etc., or a combination of the above 
cores. Each one of the local stores is a storage area 
associated with a particular SPU. In one embodiment, each 
SPU can configure its local store as a private storage 
area, a shared storage area, or an SPU may configure its 
local store as a partly private and partly shared storage. 

For example, if an SPU requires a substantial amount 
of local memory, the SPU may allocate 100% of its local 
store to private memory accessible only by that SPU. If, 
on the other hand, an SPU requires a minimal amount of 
local memory, the SPU may allocate 10% of its local store 
to private memory and the remaining 90% to shared memory. 
The shared memory is accessible by PU 4310 and by the other 
SPUs. An SPU may reserve part of its local store in order 
for the SPU to have fast, guaranteed memory access when 
performing tasks that require such fast access. The SPU 
may also reserve some of its local store as private when 
processing sensitive data, as is the case, for example, 
when the SPU is performing encryption/decryption. 



Docket No. AUS920030703US 1 49 



Atty. Ref. No. IBM-3064 



Figure 44 is a block diagram illustrating a processing 
element having a main processor and a plurality of 
secondary processors sharing a system memory. Processor 
Element (PE) 4405 includes processing unit (PU) 4410, 
which, in one embodiment, acts as the main processor and 
runs the operating system. Processing unit 4410 may be, for 
example, a Power PC core executing a Linux operating 
system. PE 4405 also includes a plurality of synergistic 
processing complex's (SPCs) such as SPCs 4445 through 4485. 
Each SPC includes a synergistic processing unit (SPU) that 
act as secondary processing units to PU 4410, a memory 
storage unit, and local storage. For example, SPC 4445 
includes SPU 4460, MMU 4455, and local storage 4459; SPC 
4465 includes SPU 4470, MMU 4475, and local storage 4479; 
and SPC 4485 includes SPU 4490, MMU 4495, and local storage 
4499. 

In one embodiment, the SPUs process data under the 
control of PU 4410. The SPUs may be, for example, digital 
signal processing cores, microprocessor cores, micro 
controller cores, etc., or a combination of the above 
cores. In one embodiment, each one of the local stores is 
a storage area associated with a particular SPU. Each SPU 
can configure its local store as a private storage area, a 
shared storage area, or an SPU's local store may be partly 
private and partly shared. 

For example, if an SPU requires a substantial amount 
of local memory, the SPU may allocate 100% of its local 
store to private memory accessible only by that SPU. If, 
on the other hand, an SPU requires a minimal amount of 
local memory, the SPU may allocate 10% of its local store 



Docket No. AUS920030703US1 



50 Atty. Ref. No. roM-3064 



to private memory and the remaining 90% to shared memory. 
The shared memory is accessible by PU 4410 and by the other 
SPUs. An SPU may reserve part of its local store in order 
for the SPU to have fast, guaranteed access to some memory 
when performing tasks that require such fast access . The 
SPU may also reserve some of its local store as private 
when processing sensitive data, as is the case, for 
example, when the SPU is performing encryption/decryption. 

The MMUs are responsible for transferring data between 
an SPU's local store and the system memory. In one 
embodiment, an MMU includes a direct memory access (DMA) 
controller configured to perform this function. 

Each SPC may be set up to perform a different task, 
and accordingly, in one embodiment, each SPC may be 
accessed using different instruction sets. If PE 4405 is 
being used in a wireless communications system, for 
example, each SPC may be responsible for separate 
processing tasks, such as modulation, chip rate processing, 
encoding, network interfacing, etc. In another embodiment, 
each SPC may have identical instruction sets and may be 
used in parallel to perform operations benefiting from 
parallel processes. 

The shared portion of the SPUs' local stores may be 
accessed by PU 4410 as well as by the other SPUs by mapping 
each shared region to system memory 4420. In one 
embodiment, PU 4410 manages the memory map for the common 
system memory 4420. The memory map table may include PU 
4410 's L2 Cache 4415, system memory 4420, as well as the 
SPUs' shared local stores. 
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A portion of system memory 4420 as shown is occupied 
by the operating system (OS 4425) System Memory 4425 also 
contains data 4440, which represents data to be processed 
by SPU 4410 as well as by the SPUs. In one embodiment, a 
process executing on the PU receives a request for a task 
involving the processing of large data. The PU first 
determines an optimtim method for performing the task as 
well as an optimum placement of the data in common system 
memory 4420. The PU may then initiate a transfer of 'the 
data to be processed from disk 4435 to system memory 4420. 
In one embodiment, the PU arranges the data in system 
memory 4425 in data blocks the size of the registers of the 
SPUs. In one embodiment, the SPUs may have 128 registers, 
each register being 128 bits long. 

The PU then searches for available SPUs and assigns 
blocks of data to any available SPUs for processing of the 
data. The SPUs can access the common system memory 
(through a DMA command, for example) transfer the data to 
the SPUs' local store, and perform the assigned operations. 
After processing the data, the SPUs may transfer the data 
(using DMA again, for example) back to common system memory 
4420. This procedure may be repeated as SPUs become 
available until all the data blocks have been processed. 

Figure 45 is a flowchart illustrating a method for 
loading data from the disk to the common system memory. 
Processing begins at 4500 whereupon, at step 4510, a task 
request is received by an executing application. The 
location of the data to be processed on disk 4540 is also 
received. The data may be a large matrix equation, for 
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example, and the requested task may be to obtain a solution 
to the matrix equation. 

At step 4515, an optimiim method for performing the 
requested task is determined. In addition, an optimiim 
block size for dividing the data is also determined. In 
one embodiment, the block size is chosen to be the size of 
the registers of the SPUs in anticipation of the parallel 
processing of the data by the SPUs. 

In step 4520, the first data block is selected, and at 
step 4525, the first data block (block 4551, for example) 
is loaded in data 4550 region in system memory 4545. A 
determination is then made as to whether more data blocks 
remain on disk 4540 requiring transfer into common system 
memory 4545 at decision 4530. If there are no more blocks 
of data to be transferred, decision 4530 branches to "no" 
branch 4534 and processing ends at 4599 . 

If there are more data blocks to be transferred, 
decision 4530 branches to "yes" branch 4532 whereupon, at 
step 4535, the next data block is loaded from disk into the 
common system memory. Processing then loops back to 
decision 4530 to determine whether there are more data 
blocks requiring transfer. 

Figure 46 is a flowchart illustrating a process for 
parallel processing data in a common system memory with a 
plurality of processors. Processing begins at 4600 
whereupon, at step 4620, the PU determines a set of SPU 
operations for performing the requested task, and at step 
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4625, the PU creates a table of a set of operations for 
completing the requested task. 

At step 4630, the PU determines an available SPU and 
sends a request to the available SPU to process a block of 
data. In one embodiment, the PU may send a request to the 
SPU by placing an appropriate value in the SPU's mailbox-a 
region of SPU memory that is continuously monitored by the 
SPU for assigned tasks. 

At step 4645, the SPU transfers the block of data to 
the SPU's local store. In one embodiment, the SPU may 
transfer the block of data using a DMA command. At step 
4650, the SPU loads the data block into the SPU's 
registers, and the SPU processes the data according to 
instructions also received from the PU. At step 4655, the 
SPU transfers the processed data block back to the common 
system memory. In one embodiment, the SPU may do so using 
a DMA command. 

A determination is then made as to whether more block 
operations are pending at decision 4660. If more block 
operations are pending, decision 4660 branches to "yes" 
branch 4662 whereupon processing loops back to step 4630 
where more SPUs are assigned data blocks for processing. 

If no more block operations are pending, decision 4660 
branches to "no" branch 4664 whereupon another 
determination is made as to whether a solution to the 
assigned task has been reached at decision 4665. If a 
solution has not yet been reached, decision 4665 branches 
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to "no" branch 4664 whereupon processing loops back to step 
4620 where a new set of SPU operations is determined. 

If a solution has been reached, decision 4665 branches 
to "yes" branch 4662 whereupon, at step 4670, the PU 
finalizes the processing. The PU may, for example, compute 
the final solution to the task by using data from all the 
processed data blocks. Processing ends at 4699. 

Figure 47 is a block diagram illustrating the 
creation, from a system of linear equations, of an 
equivalent augmented matrix. Box 4710 shows the original 
system of linear equations. A system of linear equations 
includes n unknown variables ( jc's ) linearly related by a set 
of n equations. Each variable in each equation has a 
coefficient (a's), and each equation includes a constant 
term {b's) . To solve the system of linear equations, a 
value must be found for each of the unknown variables such 
that all the equations in the system are satisfied. A 
unique solution to the system of linear equations is 
guaranteed to exist unless the determinant of the system's 
equivalent matrix (see discussion below) is zero. 

Box 4715 shows how the system of linear equations may 
be thought as an equivalent matrix equation. The matrix 
equation shown is equivalent to the system of linear 
equations since a solution to the matrix equation is also a 
solution to the system of linear equations. As shown in 
Box 4720, the matrix equation can be written in the simple 
form, ax = b, where a represents a matrix of all the 
coefficients of the unknown variables, x is a single-coliamn 
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vector of the unknown variables, and b is a single-column 
vector of the constants. 

Block 4725 shows how the coefficient matrix may be 
combined with the constant vector to yield the augmented 
matrix. In order to solve the system of linear equations, 
matrix operations are applied to the matrix equation such 
as replacing rows and columns by linear combinations of 
other rows and columns. To keep the resulting matrix 
equation equivalent (having the same solution as the 
previous matrix equation) , the same matrix operations 
should be applied to the constant vector, b . Thus, it is 
more convenient to apply these matrix operations to the 
augmented matrix, such that the operations are also applied 
to the constant vector, b. 

Figure 48 is a block diagram illustrating division of 
the linear equations coefficients into data blocks and the 
loading of the data blocks into a common memory. In one 
embodiment, the data blocks are equal in size to the 
registers of the secondary processors (SPUs) . Box 4810 
shows the undivided augmented matrix, and box 4815 shows 
the augmented matrix divided into coefficient blocks. In 
the example shown, each one of the blocks contains four 
coefficients each. 

The coefficient blocks are then loaded in common RAM 
4860. Common RAM 4860 may be accessed by the one or more 
processors of the system, which facilitates the sharing of 
data among the one or more processors. In one embodiment, 
coefficients from the same block are loaded into 
neighboring positions in common RAM 4860 in order to 
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facilitate the processing the matrix coefficients on a 
block-by-block basis in a multi-processor environment. For 
example, the coefficients from block 4820 are loaded into 
memory range 4840, the coefficients from block 4825 are 
loaded into memory range 4845, the coefficients from block 
4830 are loaded into memory range 4850, the coefficients 
from block 4835 are loaded into memory range 4855, etc. 

Figure 49 is a table illustrating examples of matrix 
operations that may be used to solve the matrix equation 
and thus the system of linear equations. The matrix 
operations in table 4900 may be created by one of the 
processors (such as a processor designated as the main 
processor) using one of the methods for solving a matrix 
equation. 

In one embodiment, the matrix equation may be solved 
by LU decomposition. LU decomposition involves factoring 
the coefficient matrix. A, into the product L-U where L 
is a lower diagonal matrix and U is an upper triangular 
matrix. A solution may then be easily obtained by solving 
for the vector U-x in the equation L-(U-x) = b and then 
solving the V-x equation for x. The above method is also 
called backward- forward substitution. The LU decomposition 
algorithm is well-known. 

Another method for solving a system of linear 
equations is Gauss elimination. The Gauss elimination 
method involves repeatedly transforming the matrix, by 
applying matrix operations, into equivalent matrices until 
the matrix is upper triangular. An upper triangular matrix 
has elements that are equal to 0 everywhere but the 
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elements along the diagonal and the elements above the 
diagonal. The unknown variables may then be easily 
computed from the upper diagonal matrix using back 
substitution. Matrix transformations typically involve 
replacing a row or column with a linear combination of the 
row or column and any other row or column. Such linear 
transformations always yield equivalent matrices— matrices 
whose solutions are the same as those of the original 
matrix. 

Column 1 of table 4900 contains a list of such 
transformations/operations. Column 2 contains, for each of 
the matrix operations, a list of blocks that contain rows 
or columns that are affected by the operation. Colxamn 3 
contains a record of whether the operations have been 
applied to the particular block, and Column 4 contains a 
list of which SPU is processing or has processed the 
particular block/operation. Free SPUs returning for a new 
assignment can be reassigned using the information in 
Column 4 to determine pending operations for a particular 
block. 

Figure 50 is a block diagram illustrating how the SPUs 
access the common memory to perform matrix operations on 
blocks of coefficients. Figure 50 shows a system having a 
main processor, PU 5010, and a number of secondary 
processors, such as SPU 5010, SPU 5015, SPU 5020,..., and SPU 
5025. For example, SPU 5010 may be accessing block 5050, 
SPU 5015 may be accessing block 5040, SPU 5020 may be 
accessing block 5055, and SPU 5025 may be accessing block 
5045. 
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Figure 51 is a flowchart illustrating the receiving of 
the linear equations coefficients and the loading of the 
coefficients into the common memory. Processing begins at 
5100 whereupon at step 5110, the number of unknown 
variables (which is also the number equations) for a system 
of linear equations is received. At step 5115, the 
coefficients of the unknown variables for the system of 
linear equations are received, and at step 5120, the 
coefficients are arranged into a matrix form. Generally, a 
system of n linear equations can be written as: 

a„aC, + 0,2^2 + ••• + flin^n = ^ 
«21^1 + «22^2 + + ^2n^n = ^2 

: ' ' 

• • • • 

where Xi-x„ are the n unknown variables, the a's are the 
coefficients of the unknown variables, and ^>'s are the 
constant terms in each equation. The solutions to the 
above linear equations are also the solutions to the 
equivalent matrix equation, ax=b, where 
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The matrix equation is equivalent to the system of linear 
equations since the solutions to the matrix equations are 
the same as the solutions to the system of linear 
equations. By operating on the a matrix to obtain 
solutions to the matrix equation, solutions to the system 
of linear equations are also obtained. 

At step 5120, an augmented matrix consisting of the 
coefficients, a, and the constants terms, b, is formed: 
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flu ^12 ••• ^1 

"21 ^22 ^2n h 

For the matrix transformations to continue yielding 
equivalent matrices, the matrix transformations should be 
applied to the coefficients as well as the constant terms. 
Thus, the augmented matrix is better form to work with. 

At step 5130, an optimum size for the coefficient 
blocks is determined. Prior to applying matrix operations, 
the matrix is divided into a number of blocks to facilitate 
applying the matrix operations by multiple processors. The 
size of the block may depend on the size of the matrix, the 
method chosen to solve the matrix, the number of available 
SPUs, etc. The size is chosen to yield the most efficient 
solving of the matrix operation. In another embodiment, 
the size of the data block may be chosen to be the size of 
the registers of the SPUs to facilitate the processing of 
the data blocks by the SPUs. 

At step 5135, the matrix is divided into blocks 
according to the determination made at step 5130, and at 
step 5140, the coefficient blocks are loaded into a common 
memory. In one embodiment, the coefficients blocks are 
loaded sequentially. The common memory is accessible by 
all the SPUs that will be sharing the task of solving the 
matrix equation. Processing ends at 5199. 

Figure 52 is a flowchart illustrating the PU 
determining a set of matrix operations to solve the linear 
differential equations. In one embodiment, the PU may use 



Docket No. AUS920030703US 1 60 Atty . Ref . No. IBM-3064 

the Gauss elimination method to obtain a solution to the 
matrix equation. The Gauss elimination method involves 
transforming the matrix into upper triangular form by 
replacing a row or a column by a linear combination of the 
row or column and one or more of the other rows or columns 
respectively. At each stage, an equivalent matrix is 
formed: one whose solution is the same as the one for the 
previous matrix. 

Processing begins at 5200 whereupon, at step 5210, the 
PU analyzes the coefficients stored in the common memory, 
and at step 5215, the PU determines the optimum method for 
solving the linear equations. For example, the PU may 
determine that Gauss elimination is the best method. 

At step 5220, the PU determines a set of matrix 
operations for solving the matrix equation. At step 5225, 
the PU creates a table containing a list of all the 
determined matrix operations. The table may be created in 
order to keep track of which operations have been^ 
completed, for which block, and by which SPU. An example 
of such a table is shown in Figure 49 . 

At step 5230, the PU programs the SPUs to perform the 
matrix operations. In one embodiment, the SPUs are 
flexible processors that can be optimized for performing 
certain tasks such as applying matrix operations to 
coefficient blocks. At step 5235, the PU instructs the 
SPUs to perform the matrix operations on a block-by-block 
basis, and at step 5240, the SPUs begin performing the 
assigned tasks. More details on the processing that takes 
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place at step 5240 are provided in the flowchart 
illustrated in Figure 53. 

After the first set of matrix operations is performed, 
a determination is made as to whether a matrix solution has 
been reached at decision 5245. If a matrix solution has 
not been reached, decision 5245 branches to "no" branch 
5255 whereupon processing returns to step 5220 where the PU 
determines a new set of matrix operations. If a matrix 
solution has been reached, decision 5245 branches to "yes" 
branch 5250 whereupon processing continues to step 5260, 

At step 5260, the PU computes the solutions to the 
matrix equation (which are also the solutions to the system 
of linear equations) from the resulting matrix. Processing 
ends at 5299. 

Figure 53 is a flowchart illustrating the SPUs 
performing the matrix operations on a block-by-block basis. 
Processing begins at 5300 whereupon, at step 5310, a free 
SPU— an SPU that is not currently involved in any other 
task— reports ready to perform pending matrix operations on 
blocks of coefficients. A determination is then made as to 
whether more block operations are pending at decision 5315. 
If there are no more block operations pending, decision 
5315 branches to ''no" branch 5325 whereupon, at step 5355, 
the PU is informed that all pending matrix operations have 
been completed. Processing ends at 5399. As shown in the 
flowchart of Figure 52, if at this time, a matrix solution 
has not been obtained, the PU may generate additional 
matrix operations . 
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If there are more block operations pending, decision 
5315 branches to ''yes" branch 5320 whereupon the SPU 
identifies a pending matrix operation and indicates that 
SPU is in the process of completing the block operation. 
In one embodiment, the SPU may identify a pending matrix 
operation and indicate the operation is being performed by 
using the table of tasks shown in Figure 49. 

At step 5335, the SPU accesses the common memory and 
loads one or more of the coefficients in its assigned block 
to begin the processing. In one embodiment, one or more of 
the SPUs may access the memory through a direct memory 
access unit. At step 5340, the matrix operation is applied 
to one or more loaded coefficients. At step 5345, the 
result of the matrix operation on the one or more 
coefficients is loaded back into the common memory. By 
doing so, the result is now accessible by the PU as well as 
by the other SPUs and there is no need to transmit the 
result to the PU or to the other SPUs. 

A determination is then made as to whether there are 
more coefficients requiring processing at decision 5360. 
If there are more coefficients requiring processing, 
decision 5360 branches to ''yes" branch 5365 whereupon 
processing return to step 5335 where one or more 
coefficients are loaded from the common memory. If there 
are no more coefficients requiring processing, decision 
5360 branches to ''no" branch 5370 whereupon processing 
returns to step 5310. At step 5310, the SPU reports ready 
to perform additional sets of block operations. 
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While particular embodiments of the present invention 
have been shown and described, it will be obvious to those 
skilled in the art that, based upon the teachings herein, 
changes and modifications may be made without departing 
from this invention and its broader aspects and, therefore, 
the appended claims are to encompass within their scope all 
such changes and modifications as are within the true 
spirit and scope of this invention. Furthermore, it is to 
be understood that the invention is solely defined by the 
appended claims. It will be understood by those with skill 
in the art that if a specific number of an introduced claim 
element is intended, such intent will be explicitly recited 
in the claim, and in the absence of such recitation no^ such 
limitation is present. For a non-limiting example, as an 
aid to understanding, the following appended claims contain 
usage of the introductory phrases ''at least one" and ''one 
or more" to introduce claim elements. However, the use of 
such phrases should not be construed to imply that the 
introduction of a claim element by the indefinite articles 
''a" or ''an" limits any particular claim containing such 
introduced claim element to inventions containing only one 
such element, even when the same claim includes the 
introductory phrases "one or more" or ''at least one" and 
indefinite articles such as "a" or "an"; the same holds 
true for the use in the claims of definite articles. 



